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We report a measurement of CP violating asymmetries in B^{B^) — s- K'^ K~ K% decays with a 
time-dependent Dalitz approach. This analysis is based on a data sample of 657 x 10^ BB pairs 
accumulated at the T(45') resonance with the Belle detector at the KEKB asymmetric-energy e+e~ 
collider. As the result of an unbinned maximum likelihood fit to the selected candidates, the mixing- 
induced and direct CP violation parameters, and Acp are obtained for — > <jf>(f020)A'g, 
-B" /o(980)iS's and other K~ Kg decays. We find four solutions that describe the data. 

There are 

(j>f{B° ^ (I}(1Q2Q)K's) = (32.2±9.0±2.6±f.4)°; 
cj>f[B"'^(l){1020)Ks) = (26.2±8.8±2.7±f.2)°; 
<j>f{B^' <}){IQ2Q)K%) = (27.3 ± 8.6 ± 2.8 ±f. 3)° and 
(/.f (B" ^ (^(1020)ifs) = (24.3 ±8.0 ±2.9 ±5.2)°. 

The values for the CP violating phase in B''^ 0(lO2O)-ft's are similar but other properties of 
the Dalitz plot are quite different for the four solutions. These four solutions have consistent <j!>i* 
values for all three B meson decay channels and none of them deviates significantly from the values 
measured in _B — >■ {cc)K'^ decays with the currently available statistics. In addition, we find no 
significant direct CP violation. 

PACS numbers: 12.15.Hh. 13.25.Hw 



CP violation in the quark sector is described in the standard model (SM) by the Kobayashi-Maskawa (KM) 
theory [l]. In this theory, the existence of a single irreducible phase gives rise to CP violating asymmetries in the 

time-dependent rates of and B decays into a common CP eigenstate Q . Specifically, for neutral B meson decays 
dominated by 6 — ccs transitions such as B'^{B ) — J/ipK^, we can measure CP violating quantities that determine 
the (f>i angle of the Unitarity Triangle. The measurements have been performed by Belle and BaBar Q 
collaborations, and provide a precise reference value for </>! because of the very small theoretical uncertainty. 

Recently, measurements of time-dependent CP violation of 6 — > s penguin-mediated B decays have become inter- 
esting because these decay modes proceed via loop diagrams and are, therefore, expected to be sensitive probes of 
the physics beyond the SM. In these decay modes, searches for new physics effects are carried out by investigating 
deviations of CP violating parameters from those determined by 6 — > ccs processes Q. 

Among these B^ decays, 5° — > K^K^Kg is one of the most promising modes because of its very small Cabibbo- 
suppressed tree diagram contribution. Previous Belle measurements of the CP violating asymmetries have been 
performed separately in the K^K~ mass region around the (/)(1020) mass and at higher K'^K~ masses jsl, while 
neglecting interference between intermediate states. It is, however, expected that the sensitivity to CP violating 
parameters would improve in a measurement using the time-dependent Dalitz plot distribution because of the correct 
treatment of interferences between various resonant and nonresonant B'^ — !■ K'^ K~ processes. 

In the decay chain T(4S)— > — > {K'^K~Kg)fti,g, where one of the B mesons decays at time t^-cc to the final 

state K^K^Kg and the other decays at time itag to a final state /tag that distinguishes between _B" and 5°, the 
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TABLE I: Summary of the contributions in the signal model. Here L is the orbital angular momentum. All fixed parameters 
are taken from Ref. [l^ except for those of the /o(980) 14] and /x T^]. 



Resonances Fixed parameters (MeV) Resonance shape L 



/o(980) M = 965±10 Flatte 

g,r=165±18 
SK=(4.21±0.33)(7^ 
(^(1020) M = 1019.455±0.020 RBW 
r = 4.26±0.04 
/x M = 1524±14 RBW 

r = 136±23 
Xco M = 3414.75±0.31 RBW 

r = 10.4±0.7 
NR no fixed parameters e~""° 
NR no fixed parameters 6"""+ 
sjNK no fixed parameters e~"*" 



{K+K- 
{K-K° 







decay rate has a time dependence given by 



|^ig(Ai,g)p = — m^ + m (1) 

4tbo L 

- gd^p- |^|2)cos(ATO<iAi) 
+ 2qIm{AA*)si-a{l\radM) , 



where r^o is the neutral B meson lifetime, Am^ is the mass difference between the two neutral B mass eigenstates, 
At ~ t^cc ~ hag, A{A) is the total amplitude of B^{B ) Kg, and the 6-flavor charge q = +1 (—1) when 

the tagged B meson is a B^{^). The Dalitz plot variables s+, s_, and sq are defined as s± = {p± +po)^, and 
So = {p+ where p+, p_, and po are the four-momenta of the K'^ , K~ , and Kg, respectively. These variables 

satisfy s+ + s- + sq = + 2to^+ + by energy-momentum conservation. In the isobar approximation [9,], the 

total amplitude for B°(B ) ^ K+K-K^ is given by the sum of the decay channels with that final state, 

where = aiC*** is a complex coefficient describing the relative magnitude and phase for the i-th decay channel, 
including the weak phase dependence. The Dalitz-dependent amplitudes, Fi{s+,S-), contain only strong dynamics 
and, thus, Fij(s+,s_) = Fi(s_,s+). The amplitudes of the contributions considered in the i?° — >■ K^K~K^ decay 
are summarized in Table HI We use the same formalism as references 13, lH ■ We utilize Flatte [12] and relativistic 
Breit Wigner (RBW) [13] hneshapes to describe the resonances. 

In the Dalitz-dependent amplitudes, A (Eq. [2]), we choose a convention in which the B'^B mixing phase {q/p) is 
absorbed into the ^ decay amplitude, a[. These complex coefficients, a'^ and a'^, can be redefined as 

in which case a resonance, i, has a direct CP violating asymmetry given by 

where the q's are restricted by definition to lie between —1 and 1. 

For cases where the contribution i is a CP eigenstate, the mixing-induced CP violating parameter, equals 
the fitted parameter di, 

0f«-^^^ = rf. (5) 
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and is related to the mixing-induced CP violating asymmetry as 

, — 2Im(a'a'*) 1 — c? „fr , , 

where rji is the CP eigenvalue of the final state. Note that Acpii) and S{i) are restricted by these definitions to lie 
in the physical region. 

This time-dependent Dalitz measurement of CP violating parameters in B'^ — > K'^K^Kg decays is based on a 
large data sample that contains 657 x 10^ BB pairs, collected with the Belle detector at the KEKB asymmetric-energy 
e+e" (3.5 on 8 GeV) collider ^16|] operating at the T(4S') resonance. The T(4S') is produced with a Lorentz boost 
factor of (3'y — 0.425 along the z-axis, which is antiparallel to the positron beam direction. Since the BB pairs are 
produced nearly at rest in the T(4S') center-of-mass system (cms), At is determined from Az, the distance between 
the two B meson decay vertices along the z-direction: At = Az/c/3j, where c is the speed of light. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon vertex detector (SVD), a 50- 
layer central drift chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), a barrel- like arrangement 
of time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) 
located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux-return located 
outside the coil is instrumented to detect mesons and to identify muons (KLM). The detector is described in 



detail elsewhere 17[. Two inner detector configurations were used. A 2.0 cm radius beam pipe and a 3-layer silicon 
vertex detector were used for the first sample of 152 x 10^ BB pairs, while a 1.5 cm radius beampipe, a 4-layer silicon 
detector and a small-cell inner drift chamber were used to record the remaining 505 x 10^ BB pairs fl8|. 

We reconstruct B candidates from an oppositely-charged kaon pair and a Kg candidate. The charged kaons are 
selected from the charged tracks having their impact parameters consistent with coming from the interaction point 
(IP). To suppress background from particle misidentification, charged tracks that are positively identified as pions, 
protons, or electrons are excluded. The particle species are identified by using particle information (PID) from the 
CDC, ACC, TOF, and ECL systems. We reconstruct Kg candidates from pairs of oppositely charged tracks having 
invariant mass within 12 MeV/c^ of the Kg mass. The direction of the Kg momentum is required to be consistent 
with the direction of vertex displacement with respect to IP ^19.] . 

We combine the K^K^ pair and Kg to form a neutral B meson. Signal candidates are identified by two kinematic 
variables defined in the cms: the beam-energy constrained mass Mbc = \/^boam ~ (~?-b)^ and the energy difference 
AE = Eb — E\^eam, where E'beam = is the cms beam energy, and b and Eb are the cms three momentum 

and energy of the reconstructed B meson candidate, respectively. We use candidates in a signal region defined as 

a 3cr ellipse around the Mbc and AE mean values: ^(^MeV/c^V^ ~'~ (45'mcV)'' where Mbo is the nominal neutral 
B meson mass [l3l|. A larger region in Mbc and AE, 5.20 GeV/c^ < Mbc and -0.30 GeV < AE < 0.50 GeV, is 
used to determine the signal and background fractions. The sideband regions used for the continuum background 
study are defined as 5.20 GeV/c^ < Mbc < 5.26 GeV/c^ and -0.10 GeV < AE < 0.50 GeV for the At distribution, 
5.24 GeV/c^ < Mbc < 5.30 GeV/c^ and -0.10 GeV < AE < 0.10 GeV excluding the rectangular region of 5.268 
GeV/c^ < Mbc < 5.30 GeV/c^ and -0.05 GeV < AE < 0.05 GeV for the Dalitz distribution. 

The dominant source of background is continuum e^e~ qq (q — u,d,s, and c) production. To reduce it, we 
require that | cos 0th | < 0.8, where 6th is the angle between the thrust axis of the B candidate and that of the rest 
of the event. This requirement retains 83% of the signal while 79% of the continuum events are removed. The BB 
background is found to be mostly originating from b c _B-decays, which peaks in the signal region with an estimated 
yield of 60 events: B^ D~[KgK~]K~^ [20| and B^ — >■ J/^pKg decays. There are also potential backgrounds from 
B° — >• D~[KgK~]K^ and — )■ 'D^[K~^K~]Kg. Backgrounds due to K-tt misidentification are also found. AU these 
peaking background decays are suppressed to a negligible level by applying ^ 1.5cr (J/'0) and ~ 2.5a (other modes) 
vetoes on the invariant masses; these vetoes are summarized in Table |TT1 For backgrounds that arise from misidentified 
particles, the invariant masses are recalculated by assuming an alternate mass hypothesis for the charged kaon. The 
remaining contribution is included in the nominal fit as the BB background component. Yields for signal, continuum 
and BB backgrounds as well as PDFs for those are described in more detail later. 

We identify the flavor of the accompanying B meson from inclusive properties of particles that are not associated with 
the reconstructed B'^ — > K^K^Kg candidate. The algorithm for flavor tagging is described in detail elsewhere [2]| . 
To represent the tagging information, we use two parameters, q defined in Eq.([T]) and r. The parameter r is an 
event- by- event Monte Carlo (MC) determined flavor-tagging quality factor that ranges from r = for no flavor 
discrimination to r = 1 for unambiguous flavor assignment. It is used only for sorting data into seven intervals. The 
wrong tag fractions for the seven r intervals, wi {I — 0,6), and the difference in w between and decays, Awi, 



are determined from data 2l|. The vertex position for the K~^K Kg decay is reconstructed using the charged kaon 
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TABLE II: Summary of the charm vetoes applied to B — >■ K^K Kg candidates. The subscript in the vetoed region indicates 
that an alternate mass hypothesis has been applied to the kaon candidates used to calculate the invariant mass term. 



Vetoed mode Vetoed region 









\M{K'l,K' 


) - Mr,- 1 


< 15 MeV/c^ 


B° - 


J/tplK+K 


-\K% 


\M{K+K- 


) - Mj/^ 


< 15 MeV/c^ 

< 15 MeV/c^ 




-^d:[k%k- 


]K+ 


\M{KlK- 




B° 




]K% 


\M{K+K~ 


)-M^o| 


< 15 MeV/c^ 




D-[KIti- 


]K+ 


\M{K%K-] 




< 15 MeV/c^ 






]K% 


\M{K+K- 


)^ - A%0 


< 15 MeV/c^ 




FIG. 1: Signal enhanced total projections of (a) AE with 5.272 GeV/c^ < M^c < 5.288 GeV/c^ (b) M^c with |A£;| < 
0.045 GeV and (c) I in the (Ai?,Afbc) signal region for the B" — >■ K'^ K~ Kg candidate events. The solid curves show the fit 
projections, the hatched areas show the continuum background component and the dotted curves show the total background 
contribution. The points with error bars are the data. 



pair and the transverse components of IP. The vertex position of /tag is obtained using tracks that are not assigned 
to the K+R-R'^ candidate and IP. 

We find that 1.5 % of the selected events have more than one B'^ — > K^K^Kg candidate. In these events, we 
choose the B'^ candidate that is formed from the most kaon-like charged kaon candidate and the Kg candidate closest 
to the nominal Kg mass. 

After all the selections are applied, we obtain 98982 candidates in the Afbc-Ai? fit region, of which 2333 are in the 
signal region. We extract the signal yield using a three-dimensional extended unbinned maximum likelihood fit to 
the distributions of AE, Afbc and the flavor-tag quality (r) interval, I, for the selected B^ — >■ K~^K~Kg events. For 
the probability density function (PDF) of the signal component, we use a sum of two Gaussians (a single Gaussian) 
for the AE (A/bc) shape. All parameters of the PDFs are free in the fit, except the ratio of the area of the broader 
Gaussian component to that of the core Gaussian, and the width of the broader Gaussian in AE. These additional 
parameters are fixed from the results of a fit to a i?'' — D~ [Kg7T~]Tr'^ data control sample. For the continuum 
background component, the AE (A/bc) shape is modeled by a first-order polynomial (an ARGUS [22!|) function, with 
shape parameters floated in the fit. The BB background component is parameterized by two-dimensional binned 
histograms from MC. In the fit, the total signal, continuum and BB background yields are also free parameters. 
The fit yields 1176 ± 51 signal events in the signal region. The projections of the AE, A/bc and I distributions for 
the candidate events are shown in Fig. [T] The average signal, continuum and BB fractions in the signal ellipse are 
calculated to be ~50 %, ~49 % and ~1 %, respectively. The event-by-event signal probabilities as a function of AE, 
A/bc and / obtained with this fit are used in the unbinned maximum likelihood fit with a time-dependent Dalitz 
approach that is used to extract CP violation parameters. 

In a Dalitz plot as a function of (s+, s_), signal and continuum events densely populate the kinematic boundaries 
with low So, which correspond to the (/)(1020) and /o(980) resonances. Large variations in a small area of the Dalitz 
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(a) (b) 




FIG. 2: The Dalitz distribution based on our signal model of GEANT-based signal MC (a) with the normal Dalitz parame- 
terization, s+ and s_, and (b) with the square Dalitz parameterization, m' and 9' . The dashed red boxes indicate the regions 
where most of the signal components and the background are located. 



plane make it difficult to use histograms to describe the background. Therefore, we apply the transformation, 

ds+ds^ \AetJ\dm' d0' , (7) 
where J is the Jacobian of this transformation. The parameters m' and 9' are given by the transformation. 



— arccos 2 



max 




Tinin 
'■0 



1 



and 



1 



(8) 
(9) 



where toq is K~ invariant mass, m}^^^ and m™'" are kinematic limits of mg, is the helicity angle, defined as 
the angle between the if ^ and the in the K^K^ rest frame. With this transformation, the Dalitz plot turns into 
a "square Dalitz plot" with a smooth density variation. Figure [2] (a) and (b) show the Dalitz distributions based on 
our signal model with the usual Dalitz parameterization, s+ and s_, the square Dalitz parameterization, m' and 9\ 
respectively. As can be seen, the highlighted region where most of the signal and background events are located, is 
magnified in the square Dalitz parameterization. The square Dalitz plot is described in detail elsewhere 13, 23 1. 
The PDF expected for the signal distribution, T'sig, is given by 



P,ig(m', 9', At, q) = e(m', (?')l^igK, 0\ At, q)\' 



R. 



■Slg7 



where 



\A,ig{m',9',At,q)\' 



\detJ\- 



\^t\/Tg„ . 



il~qAwi)i\A\^ + \An 



(10) 



(11) 



Atbo 

- qil-2wi){\Af ^\Af)cosiAmdAt) 
+ 2q{l -2wi)Im{AA*)sin{AmdAt) , 

which accounts for CP dilution from the incorrect flavor tagging. This function is convolved with the At resolution 
function i?sig [7|; the impact of detector resolution on the Dalitz plot is ignored because the intrinsic widths of the 
dominant resonances are larger than the mass resolution. We determine the variations of the signal detection efficiency 
across the Dalitz plane due to detector acceptance, e(m', 9'), by using a large MC sample. 
The PDF for continuum background is 



Vaaim' 



.^t,q) 



1 + qA,^{9') 



Hgg{m\9')Pgg{At), 



(12) 
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where Hqq, Aqq, and Pqq{At) are the Dahtz distribution PDF, the Dahtz-plot-dependent flavor asymmetry and the 
At PDF, respectively. The function Pqq is modeled as a sum of exponential and prompt components, and is convolved 
with a double Gaussian that represents the resolution. All parameters of Pqq are determined by a fit to the At 
distribution in the sideband region that is defined above. The Dalitz distribution PDF, Hqq^ is a two-dimensional 
binned histogram PDF. To determine the PDF, we use the sideband region around the signal region with a less 
restrictive requirement, | cos(?th| < 0.92, to increase statistics. We have checked that the Dalitz distribution for this 
sideband region is similar to that for the signal region, using a MC sample. There is a flavor asymmetry Aqq due to 
the jet-like topology of continuum because a high momentum K^(K~) in f^cc is accompanied by a high momentum 
K^{K^) in /tag! to account for this, we extract the Dalitz plot asymmetry using almost the same region as the region 
used in Hqq extraction: since we find no correlation between 0' and Mbc, we enlarge the lower limit of the sideband 
region in this fit from 5.24 GeV/c^ to 5.2 GeV/c^ in Mbc. 

Using high-statistics MC sample, we find no CP violating asymmetry in the background coming from charmless 
and charmed B decays. Therefore, the PDFs for B^B^ and B^B backgrounds are given by 



VB^B-{T^\e\At) = HB+B-{m\e')\AB^B-{^t)\'' ® Rb+b- , 

Vj.,-go{m',e',At) = Hg„go{m' ,9')\Ag„go{At)\'' ® Rg„^o, 



(13) 
(14) 



respectively. Dalitz distribution PDFs as well as Hg^g_^g„-gOy are modeled with two-dimensional histograms from 
MC. The At PDF for both models, I^^+^-j-^ob")!^! ^-re described by exponential functions with effective lifetimes 
while Rb+b~{b'>'b°) ^""^ ^^"^ resolution functions. The effective lifetimes are obtained from fits to the MC sample. 

To account for a small fraction of events with large At values not yet described by either signal or background 
PDFs, an outlier PDF is introduced, Vol — HoiG{At), where G is a Gaussian and Hoi is the two-dimensional binned 
histogram PDF of the Dalitz plot of data itself. 

For the j-th event, the following likelihood function is evaluated: 



P,{m',e',At,q;AE,Mkc,l) = (1 - foi) 



. k 



fkiAE,A'h„l)Vk{m',9',At,q) 



+ foiVoi{m',9',At) (15) 



where k runs over a total of four components including signal and backgrounds. The probability of each component 
{fj) is calculated using the result of the AE-M^c-l fit on an event-by-event basis. 

As there is only sensitivity to the relative amplitudes and phases between decay modes, we fix a(K+ K-)^Yt ^ 
and b(^K+K-)fiu_ — 0°- addition, /x and non-resonant contributions are combined and have a single common CP 
violating parameters. The combined component is referred to as "others" throughout this paper. The parameters 
d^^g and c^^g are fixed to the world average b — ccs values of 21.5° and 0, respectively. We determine 19 parameters 
of the Dalitz plot and CP asymmetries by maximizing the likelihood function C = Y[j Pji where the product is over 
all events. 

We find four preferred solutions with consistent CP parameters but significantly different amplitudes for /o(980)ii'5 
and fx Kg. The fitted results are summarized in Table Hill These are obtained by performing a large number of fits 
with random input parameters. For each resonance, i, the relative fractions can be calculated as 



(|a,P + |a,p) / F,{s+,s^)F*{s+,s^)ds+ds 
f{\A\^ + \W)ds+ds. 



(16) 



where the sum of fractions over all decay channels may not be 100% due to interference. Table IIVI summarizes the 
relative fractions for all solutions. 

By translating the fit results using Eqs. 2] and [SJ we determine the time-dependent CP violating parameters of 
B° ^ /o(980)ifg and B° (/)(1020)A"g decays and other B° decays with the K+R-Rl final state. Table |V] 
summarizes the CP violating parameters for all solutions. 

In Table ITVl /0(io2o)icg is similar for all four solutions but //^(gso)^^ and ff^K% significantly different. These 
four solutions are due to interference between the /o(980) and non-resonant component, and interference between 
the /x and non-resonant component, and are characterized by different relative fractions for /o(980) and /x- In 
order to distinguish these solutions with the current statistics, we use external information from B^ — >■ tt^tt 
and the property that /o(980) decays to either tt+tt" or K^K~ . We calculate the branching fraction, B{B^ 
/o(980)[7r+7r~]Ar§), based on the branching fraction of _B" Tr+Tr^ATg and the relative fraction of /o(980)Ar§ in the 
^ decay 24]. Similarly, we can calculate the branching fraction, B{B^ — /o(980)[Ar+Ar~]Ar^) from 

B(/o(980)-^7r7r) calculated to be 0.47 ± 0.10 for Solution 1 



B" 



Table |IVl The fraction, //^(gso) 



S(/o(980)- 



hB{fo{9S0)^KK)' 
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TABLE III: Time-dependent Dalitz plot fit results for tfie four solutions witli statistical errors. The phases, bi and di, and a 
are given in degrees and GeV~^c*, respectively. 



Parameter 



a/o(980) 

10(1020) 



^/o(980) 
^.^>(1020) 



C/o(980) 

c<;>(io20) 



Solution 1 
29.3 7 

'-'•'J'' -0.06 

5.2 ± 0.8 
2.03 



Solution 2 
^t:^ — 



+0.31 
-0.28 

D-0 -6.7 
+4.9 
4.3 



53.0 
0.67 
7.0 
2.53 



25.9 



-19.9 
+0.10 
-0.23 
+ 1.2 
-2.5 
+0.46 
-0.89 

20.8 ° 
40.2 til 



-16.0 
-34.5 
-32.6 
-28.0 
126.5 _g2.7 
-123.5 



+ 10.0 
-13.2 
+ 14.0 
-14.5 
+8.3 
-8.9 
+22.5 
26.6 



83.4 
108.7 
-106.0 
-36.2 
113.6 



+0 — 

-7.3 
+ 15.8 
-15.5 
+ 12.9 
-13.8 
+25.0 
-28.3 
+ 13.6 
-16.2 



Solution 3 
31.8 

0.56 trot 

15.6 tli 
2 16 +°-'^^ 

10.3 tli 

+6.4 
-5.4 
flO.O 
-13.5 
f 14.1 
-14.5 
+8.7 
-8.2 
+25.3 
-27.8 
+ 19.0 
-31.4 



Solution 4 
64.1 m 

0.71 till 
23.9 +^-^ 

2.89 



29.7 



3.1 
+0.56 
-0.45 

24.3 ttl 

+ 5.9 
-6.0 



21.7 



+ 13.4 
-11.4 



-143.5 ^7.6 



-1.8 

-7.0 

92.1 
-35.7 
113.5 
-141.2 tli 



109.3 
149.2 
35.7 
45.96 
98.8 
23.7 



+ 10. S 
-8.2 
+ 17.5 
-16.6 
+6.9 
-5.9 
+22.4 
-30.2 
+23.0 
-21.8 
+22.2 
-24.6 



0.16 



+0.16 
-0.15 



-0.02 ± 0.10 



0.10 to'lr -0-01 ± 0.11 
-0.04 ± 0.09 0.01 ± 0.10 



0.09 ± 0.07 
-0.10 ± 0.09 



•^others 


0.07 ± 0.06 


0.03 ± 0.08 


nni 

'^■'^^ -0.05 


-n 02 

'J-'J^ -0.06 


rf/o(980) 
^0(1020) 
^^others 


31.3 

32.2 in 

24.9 til 


26 1 +'^-" 

^O-J^ -6.6 

26.2 til 
29.8 til 


25.6 l'^;« 

07 +8.6 
^ ' -8.0 

26.2 til 


26.3 til 
24.3 tf° 
23.8 til 


a 


12 +"-"^ 

U-J^^ -0.04 


0.06 ± 0.04 


0.10 ± 0.04 


0.18 ± 0.03 


-21og£ 


10201.7 


10198.6 


10204.5 


10208.9 



TABLE IV: Summary of the relative fractions (%), the errors are statistical only. 



Parameter 



//o (980) Kg 
f<p{1020)K^ 

./'(K+if-)NRXO 

^ )nr 

K+ 



Solution 1 Solution 2 Solution 3 Solution 4 



26.0 ± 7.4 
14.2 ± 1.2 
5.10 ± 1.39 



54.0 ± 9.6 
14.5 ± 1.2 
5.89 ± 1.86 



3.73 ± 0.74 3.71 ± 0.73 



26.4 ± 7.8 
14.2 ± 1.2 
39.6 ± 2.6 
3.68 ± 0.73 



68.1 ± 12.3 
14.4 ± 1.2 
59.0 ± 3.0 

4.15 ± 0.79 



138.4 ± 44.8 175.0 ± 52.6 157.4 ± 29.5 48.1 ± 11.7 
1.65 ± 4.17 21.0 ± 17.3 4.63 ± 6.76 7.87 ± 4.78 
26.0 ± 12.9 78.0 ± 36.2 38.6 ± 18.1 6.27 ± 3.81 



215.2 ± 47.5 352.0 ± 66.8 284.5 ± 36.3 207.9 ± 18.4 



and 0.30 ± 0.07 for Solution 2. The value of //o(98o)->7r7r is also determined by the BES Collaboration, which uses 
the same parameterization for the /o(980), to be 0.75 ± 0.12 [14,]. Therefore, the solutions with a low fo{980)Kg 
fraction (Solution 1 and 3) are preferred. It is likely that the /x described in this analysis and the i?° — > tt^tt~ Kg 
analysis is the same state, /o(1500). If this is the case, the ratio, B(fo(i50o)l^Kl<) ' calculated to be 1.0 ± 0.7 
for Solution 1 and 0.13 ± 0.09 for Solution 3. As the world average of this ratio is 4.1 ± 2.5 the solutions with a 
low fxKg fraction (Solution 1 and 2) are preferred. Table IVll summarizes these values for each solution. Altogether, 
we conclude that Solution 1 is preferred from all currently available external measurements. The mass projections 
onto (a) M{K^K+), (b) M{K^K~) and (c) M(K+K-) distributions for Solution 1 are shown in Fig. H and (a) At 
distribution and (b) raw asymmetry in the 0(lO2O)ifg region are shown in Fig. 21 The full correlation matrix is given 
is Tables Ivnl and I Villi Likelihood scans of (j)'^ for all four solutions are obtained by fixing cf)'^ and redoing the fit. 
We also perform scans that include the systematic and model errors by convolving the likelihood with a Gaussian 
with the width set to the quadratic sum of the systematic and mode uncertainties. These are shown in Figs. [5l [H 

The sources of systematic uncertainties and their contributions are summarized in Table IIXI The systematic 
errors in the vertex reconstruction include uncertainties in the IP constraint, charged track selection based on track 
helix errors, vertex reconstruction quality. At requirement, tracking error corrections, Az bias, and imperfect SVD 
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TABLE V: Time-dependent CP violating parameters for the four solutions, where the first error is statistical, the second is 
systematic and the third is the Dalitz plot model uncertainty. 





Solution 1 




Solution 2 




Solution 3 




Solution 4 






-0.30 ± 


0.29 ±0.11 


±0.09 


-0.20 ± 


0.15 ±0.( 


38 


±0.05 


±0.02 ± 


0.21 ±0.09 


±0.09 


-0.18 ±0.14 ±0.08 


±0.06 


-^f (/o(980)7^§) 


(31.3 ± 


9.0 ± 3.4 ± 


4.0)° 


(26.1 ± 


7.0 ±2.4 


± 


2.5)° 


(25.6 ± 


7.6 ± 2.9 ± 


0.8)° 


(26.3 ± 5.7 ± 2.4 ± 


5.8)° 




+0.04 ± 


0.20 ±0.10 


±0.02 


±0.08 ± 


0.18 ±0. 


10 


±0.03 


-0.01 ± 


0.20 ±0.11 


±0.02 


±0.21 ±0.18 ±0.11 


±0.05 


(0(lO2O)7^g) 


(32.2 ± 


9.0 ± 2.6 ± 


1.4)° 


(26.2 ± 


8.8 ±2.7 


± 


1.2)° 


(27.3 ± 


8.6 ± 2.8 ± 


1.3)° 


(24.3 ± 8.0 ± 2.9 ± 


5.2)° 


Acp (others) 


-0.14 ± 


0.11 ±0.08 


±0.03 


-0.06 ± 


0.15 ±0.( 


38 


±0.04 


-0.03 ± 


0.09 ± 0.08 


±0.03 


±0.04 ±0.11 ±0.08 


±0.02 


<j}f (others) 


(24.9 ± 


6.4 ± 2.1 ± 


2.5)° 


(29.8 ± 


6.6 ±2.1 


± 


1.1)° 


(26.2 ± 


5.9 ± 2.3 ± 


1.5)° 


(23.8 ± 5.5 ± 1.9 ± 


6.4)° 



TABLE VL Comparison of external information with each of the four solutions. 





//o(980)^7r7r 


B(/oC1500)-+7r7r) 
S(/o(1500)^KA-) 


Solution 1 


0.47 ±0.10 


1.0 ±0.7 


Solution 2 


0.30 ± 0.07 


0.91 ±0.64 


Solution 3 


0.46 ±0.10 


0.13 ±0.09 


Solution 4 


0.25 ±0.06 


0.09 ±0.06 


External information 


0.75 ±0.12 [14] 4.1 ±2.5 fl3] 



alignment. The parameters for flavor tagging and resolution function, physics parameters, background At shape and 
signal probability are varied by ±lcr. For each histogram, systematic errors are estimated using 100 sets of pseudo- 
experiments, generated by statistically fluctuating samples to create the histogram. Samples of pseudo-experiments 
showed some fitting bias for CP parameters due to low statistics in each sample. We take this bias as a systematic 
uncertainty. The effect of misreconstruction is accounted for by comparing the fitted results of signal MC samples with 
and without misreconstructed events. The efficiency histogram also includes systematic uncertainties in the Dalitz- 
dependent correction factors for Kg, PID and tracking efficiency. For tag-side interference [25| . pseudo-experiments 
are generated with and without tag-side interference and the difference is taken as a systematic error. The fixed 
masses and widths of the resonance form factors in the signal model are varied by their errors, the changes in the 
results are taken into account in the systematic errors. To take into account a systematic uncertainty due to mass 
resolution in the (/)(1020) mass region, the width of the 0(1020) is varied from 4.26 MeV/c^ to 5.4 MeV/c^ and the 
difference in the fitted result is taken. The systematic uncertainty due to Blatt-Weisskopf barrier factors [2^ in 
RBW is determined by taking the difference in the fitted results with and without these factors. The shape of the 
non-resonant component is empirically chosen, so different parameterizations are possible. This includes modeling 
the non-resonant part with the tail of a Breit-Wigner, i?jvfl(s; a) = ia/{s ± ia) and a power law whose exponent is a 
fitted parameter, Rtqii{s\ a) = s^". A possible variation in the model of the /o(980) Flatte function, is also considered 



TABLE VIL Statistical correlation matrix for Solutionl. 



a/o(980) 

'^i^>(1020)kO 

«{X«K+)nr 

^/o(980) 
^<A(1020)Kg 



1.00 



0.13 
1.00 









«(KOif-)NR 


&/o(980) 


^.^>(1020)if^l 






0.42 


0.13 


0.29 


0.17 


-0.04 


0.22 


-0.18 


0.14 


0.57 


0.71 


0.54 


0.71 


-0.88 


-0.35 


-0.38 


0.15 


1.00 


0.45 


0.50 


0.51 


-0.44 


-0.09 


-0.09 


0.12 




1.00 


0.41 


0.52 


-0.71 


-0.28 


-0.31 


0.37 






1.00 


0.67 


-0.51 


-0.05 


-0.12 


0.17 








1.00 


-0.73 


-0.14 


-0.29 


-0.01 










1.00 


0.38 


0.40 


-0.09 












1.00 


0.17 


-0.01 














1.00 


-0.01 
















1.00 
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FIG. 3: The mass projections onto (a) M{KsK+), (b) M{KgK~) and (c) M{K+K~) distributions (the inset shows the 
projection near the 0(1020) resonance) for the K'^ K~ K% candidate events in the signal region, using Solution 1. In 

(a-c), the red solid curves show the fit projections while the blue hatched areas and the green dashed curves show the qq and 
total background components, respectively. The points with error bars are the data. 



using a different parameterization [27[. We also include a possible contribution from the spin 2 /2(1270) resonance in 
tfie signal model. The differences in the fit results from these alternate Dalitz plot parameterizations were summed in 
quadrature. The total systematic uncertainty is obtained by summing all of the above contributions in quadrature. 

In summary, for the first time in Belle we perform a measurement of the CP violating asymmetries in B^{B ) 
K'^K^Kg decays with the time-dependent Dalitz approach. There are four solutions that describe the data well. 
These give similar values for the CP violating phase in the 6 — > s penguin mode, i?'^ — !• (j){1020)Kg, 

(l)f{B° ^ <j){1020)K°s) = (32.2 ±9.0 ±2.6 ±1.4)°; 

(jjf(^B° <j){1020)K°s) = (26.2 ±8.8 ±2.7 ±1.2)°; 

(l)f{B° ^ (j){1020)K"s) = (27.3 ± 8.6 ± 2.8 ± 1.3)° and 

(l)f{B° ^ (j){1020)K"s) = (24.3 ±8.0 ±2.9 ±5.2)°. 



11 




FIG. 4: (a) At distribution and (b) raw asymmetry for the — > K'^ K K% candidates in the 0(lO2O)A'g region, |M^+j^- — 
M^(i02o) I < 0.01 GeV/c^, with good tags, r > 0.5, using Solution 1. In (a), the blue solid and red dashed curves show the 
fitted results with i?° and ^ tags, respectively. The dotted black curve shows the background component with and 
tags. In (b), the solid curve shows the fit projection and the dashed curve shows the Standard Model expectation from the 
time-dependent CP asymmetry measurement in fe — >■ ccs decays. 



TABLE VIII: Statistical correlation matrix for Solutionl. 













f^othcrs 


C/o (980)^0 




Cot hers 


a 


a/o(980) 


0.02 


-0.09 


0.11 


0.01 


0.14 


-0.20 


-0.05 


-0.01 


-0.35 


'*<;>(io20)ko 


0.15 


-0.24 


0.22 


0.07 


0.07 


0.02 


0.01 


0.07 


-0.72 




0.11 


-0.28 


0.15 


0.08 


0.06 


0.01 


0.01 


0.05 


-0.65 




0.17 


-0.12 


0.19 


0.06 


0.08 


0.00 


0.02 


0.05 


-0.54 


«{KOif+)NR 


0.26 


-0.25 


0.20 


0.08 


0.26 


-0.14 


-0.01 


-0.05 


-0.81 




-0.16 


-0.50 


0.22 


0.07 


0.19 


-0.09 


-0.01 


-0.05 


-0.77 


^/o(980) 


-0.21 


0.29 


-0.30 


-0.11 


-0.03 


-0.09 


-0.03 


-0.16 


0.72 


^<^(1020)ifg 


-0.20 


-0.01 


-0.00 


-0.09 


0.11 


-0.20 


-0.03 


-0.09 


0.15 




-0.09 


0.01 


-0.10 


-0.04 


-0.02 


-0.00 


-0.01 


-0.03 


0.21 




0.00 


0.18 


0.11 


0.05 


0.17 


-0.07 


-0.01 


0.01 


-0.07 




1.00 


0.03 


0.01 


-0.01 


-0.05 


0.02 


0.02 


0.05 


-0.30 


^(/ft)/f-)NR 




1.00 


-0.08 


-0.02 


-0.00 


0.01 


0.02 


0.02 


0.63 


'^/o (980) Kg 






1.00 


0.60 


0.59 


0.14 


-0.02 


0.47 


-0.20 


'^<I>(1020)K° 








1.00 


0.35 


0.45 


0.08 


0.37 


-0.07 


^others 










1.00 


-0.39 


-0.09 


-0.13 


-0.14 


C/o(980)itO 












1.00 


0.24 


0.60 


0.07 


C,^>(1020)ifg 














1.00 


0.14 


0.00 


Cothcrs 
















1.00 


0.00 


a 


















1.00 



These CP violating parameters are consistent at the current level of precision with the measurement of the CP 
violating phase in & — >■ ccs processes such as B'-* — > J/ipK'^, which is (22 ± 1)°. Previous measurements used K'^K^ 
selections around the 0(1020) mass and in the higher K^K~ mass region. Here we establish a superior analysis 
procedure for obtaining CP violating asymmetries without the uncertainty from interference among different resonant 
contributions, and therefore this represents an important step toward measurements with higher statistics such as in 
Super B-factory experiments. 
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FIG. 5: Likelihood scans of (j>f {fo{980)K^) for (a) Solution 1, (b) Solution 2, (c) Solution 3, and (d) Solution 4. The solid 
(dashed) curve contains the total (statistical) error and the dotted box indicates the parameter range corresponding to ilcr. 
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FIG. 6: Likelihood scans of (j>f {(l>{W20)K°s) for (a) Solution 1, (b) Solution 2, (c) Solution 3, and (d) Solution 4. The solid 
(dashed) curve contains the total (statistical) error and the dotted box indicates the parameter range corresponding to ilcr. 
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